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Abstract 

Supernova neutrinos could be well-suited for probing neutrino decay, since decay 
may be observed even for very small decay rates or coupling constants. We will 
introduce an effective operator framework for the combined description of neutrino 
decay and neutrino oscillations for supernova neutrinos, which can especially take 
into account two properties: One is the radially symmetric neutrino flux, allowing 
a decay product to be re-directed towards the observer even if the parent neutrino 
had a different original direction of propagation. The other is decoherence because 
of the long baselines for coherently produced neutrinos. We will demonstrate how 
to use this effective theory to calculate the time-dependent fluxes at the detector. In 
addition, we will show the implications of a Majoron-like decay model. As a result, 
we will demonstrate that for certain parameter values one may observe some effects 
which could also mimic signals similar to the ones expected from supernova models, 
making it in general harder to separate neutrino and supernova properties. 
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1 Introduction 



Neutrino decay [1,2] lias been considered as an alternative to neutrino os- 
cillations, especially for atmospheric [3-8] and solar [3,9,10] neutrinos. Fur- 
thermore, sequential combinations of neutrino decay and neutrino oscillations 
have been studied {e.g., MSW-mediated solar neutrino decay [11,10]), as well 
as a combined treatment [12]. Currently, one of the most favorable models 
for neutrino decay is to introduce an effective decay Lagrangian which couples 
the neutrino fields to a massless boson carrying lepton number, i.e., a complex 
scalar field or a Majoron field [13-16]. One possibility for such a Lagrangian 
for the case of Majoron decay is 



where the z/'s are Majorana mass eigenfields and J is a Majoron field. An 
interaction Lagrangian like this in general also implies coupling to active neu- 
trino mass eigenstates. It was shown in Ref. [12] that interference effects, such 
as neutrino oscillations, are in principle possible before and after decay. There 
are several constraints to Majoron models, such as the exclusion of pure triplet 
models by the width from LEP at CERN. However, fine-tuning or small 
triplet admixtures to singlet Majorons, as well as more sophisticated models 
have been proposed to circumvent this problem. For a more detailed discussion 
see, for example, Refs. [3-5]. Alternative neutrino decay models, such as the 
one mentioned in Ref. [17], proposing decay into more than one neutrino, have 
also been suggested. Nevertheless, since we are interested in fast neutrino de- 
cay possibilities, we will focus on models for decay of one neutrino into exactly 
one neutrino, which can be described by a Lagrangian similar to Eq. (1). 

In this paper, we extend the effective operator formalism introduced in Ref. 
[12] to allow a time- dependent treatment of a supernova neutrino flux modified 
by neutrino decay and neutrino oscillations. In earlier works [18-21], different 
aspects of supernova decay have been considered, which we want to combine 
into a general formalism. For the sake of the analytical accessibility, however, 
we have to make certain simplifications. Since we want to focus on the neu- 
trino propagation in vacuum, we thus ignore matter effects as well as any 
other effect within the supernova. Nevertheless, as we will see, taking a su- 
perposition of mass eigenstates obtained from a numerical calculation within 
a supernova instead of an initial flavor state, can be done in a quite straight- 
forward way in our formalism. Furthermore, we restrict ourselves to maximal 
one intermediate decay between production and detection[^ or almost stable 

^ Repeated decays correspond to higher order processes suppressed by the small 
coupling constants. However, depending on the decay model, for the very large 
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decay products^. We are especially interested in possible modifications of the 
time dependence of a supernova neutrino signal. Therefore, we neglect other 
effects such as (cosmological and non-cosmological) redshiftQ, asymmetries 
of the supernova flux, and spatial extension of the neutrino production re- 
gionQ. However, we take into account dispersion by different traveling times 
of different mass eigenstates, since we want to study the physical effects of 
combining neutrino decay and neutrino oscillations. In addition, we incor- 
porate that neutrino oscillations are affected by wave packet decoherence at 
long traveling distances if the neutrinos are coherently produced or leaving 
the supernova as superpositions of mass eigenstates at all. Since coherence 
of supernova neutrinos is still a question under discussion, we will a priori 
assume that supernova neutrinos are produced coherently as flavor states, but 
demonstrate how to obtain the incoherent limit thereafter. When taking into 
account the MSW effect, coherent emission can be achieved for appropriate 
parameter values violating the adiabaticity condition. In addition, for a three 
flavor neutrino scenario obtaining exactly one (pure) mass eigenstate emerg- 
ing from the supernova is no longer valid. It may even be a superposition of 
two mass eigenstates, such that neutrino oscillations of supernova neutrinos 
(close to the supernova) cannot be excluded in general. Nevertheless, we will 
use the incoherent limit in most of our examples, since this is the most often 
assumed case. 

We will cover three main topics in this paper: First, in Sec. 2 we will discuss 
the loss of coherence because of the long baselines and how this is modified 
by intermediate decays. Especially, we will use Majoron decay as an example 
for demonstration. Second, in Sec. 3, we will introduce a general formalism in 
order to be able to calculate time- dependent fluxes at the detector. It will also 
implement the coherence issues mentioned before. Since it will also work for 
decay models similar to Majoron decay, it will be treated somewhat decoupled 
from the initial discussion. Third, in the remaining sections, we will give some 
applications. In Sec. 4, we will demonstrate the time smearing of a source pulse 



supernova distances they may be relevant making the geometrical discussion much 
more complicated. For the implementation of repeated decays see also Ref. [12]. 
^ Since we assume that particles other than the neutrinos, such as Majorons, are 
not detectable, the terms decay products and secondary neutrinos will refer to the 
neutrinos produced by decay. Note that in this paper the secondary neutrinos can 
be either active or sterile neutrinos. 

^ We ignore redshift, since we know from supernova statistics that the observable 
supernovas are quite close to the Earth. In addition, we assume the supernovas not 
to be moving too fast within the cosmological comoving frame. 

Compared with the distance to the Earth, supernovas can for our purpose be 
treated as point sources. In addition, oscillation phases of neutrinos produced by 
a core-collapse supernova are not necessarily averaged out by the extension of the 
production region [0(10 km)], which is typically much smaller than the oscillation 
length in vacuum [0(1000 km)]. 
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by neutrino decay. Coming back to Majoron decay, in Sec. 5, we will take a 
closer look on the dynamics of a Majoron decay model and its implications. 
In the next section, Sec. 6, we will calculate the flux at the detector by also 
taking into account path lengths of different traveling paths. In addition, we 
will apply the results for Majoron decay we obtained in Sec. 5 in order to 
show the implications for this decay model. Finally, in Sec. 7, we will show 
the peculiarity of observing interference effects at the detector in the limit of 
quite large decay rates, even if the coherence lengths of the problem are much 
shorter than the baseline length. In almost all sections, we will use simple 
examples to visualize the indicated effects. 



2 Coherence in neutrino oscillations and intermediate decays 

The coherence length for neutrino oscillations is of the order L'^°^ ~ \aS\ ' 
where is the width of a (Gaussian) wave packet [22-25] or the wave packet 
overlap in the source or detector [26]. Thus, for typical values for supernova 
neutrinos coherence seems to be destroyed on their way to the Earth. However, 
since the decay may happen close to the source or detector, the distance 
between decay and production or detection may be much shorter than the 
coherence length. We will show that we may thus see coherent interference 
effects in certain cases, e.g., for very early decays. For this, we need to take 
into account that the wave packet width, which enters the coherence length, 
is in fact = (cr^)^ + (c^)^, i-^-, the squared sum of the widths of the 
production P and detection D processes [23,26]. One can visualize this by a 
detection process working on a much longer timescale than the production 
process. Therefore, though the wave packets may be well separated to a third 
observer at arrival, the detector will not be able to resolve them due to the 
time resolution of its detection process. However, neutrino oscillations may be 
washed out, e.g., for too large > L°^'^, since the detector then averages 
over the oscillation length. For supernova neutrinos, the width is because of 
< 10~^^ m ^ (7^ determined by the detection process only. Taking into 
account intermediate neutrino decays, the detection process D corresponds to 
the decay process X with its width . For Majoron decay of a relativistic 
mass eigenstate i of energy E into a mass eigenstate j, determined by the 
Majoron coupling constant gij, one can estimate the wave packet overlap by 
the decay rate Fj^ as 



where the exact value of Tij depends on the interaction Lagrangian [27]. Thus, 
the weaker the coupling constants g^j are, the slower is the process and the 
longer is the spatial extension of the wave packet. Especially, for g^j — > the 
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wave packet width would become infinite. Nevertheless, we will see that neu- 
trino oscillations are not possible in this limit because of another constraint to 
be satisfied. Equation (2) leads to a coherence length for supernova neutrinos 



coh,X 

'ab = 



= O 



mimjAml,,gfj 



(3) 



Here Am^;, is the mass squared difference of the neutrino oscillation considered 
(before or after decay). For instance, an incoming superposition of two mass 
eigenstates, oscillating in the flavor basis, say 1^2 and 1/3, may decay into an 
outgoing mass eigenstate, say z/i, by non-zero Majoron coupling constants g^i 
and (721- Then the index i of the parent neutrino is 3 or 2, the index j of 
the decay product is 1, and the Am^^ considered is equal to Am|2- From the 
last equation, Eq. (3), we see that especially for small coupling constants the 
coherence length can be quite long. However, in Refs. [23,26], the upper bound 



was found for the coherence length, determined by the condition < L"^^ 
necessary for neutrino oscillations not to be washed out by the spatial wave 
packet extension. Thus, this implies that L'^°^ < (10^ ~ 10^^) m for Am^ ~ 
(10-^ ~ l)eV^ which is, in principle, much shorter than the typical distance 
of a supernova. In addition, in order to observe the oscillations before or after 
decay, this means for Majoron decay that 



This is a finite lifetime condition similar to the one in Ref. [24] for the decaying 
muons used for neutrino production. 

Since we want to obtain a treatment which is independent of the coherence 
discussion, we use the wave packet approach in Ref. [23] by assuming that 
coherence is destroyed by a factor exp[— (//L'^"^)^]. Here L^"^ is given by Eq. (3) 
with the respective widths of the processes considered. In addition, model- 
dependent constraints may have to be applied to L^°^. The following two 
cases will be discussed in this paper: 

(1) Coherence is lost before decay and detection, i.e., all coherence lengths 
involved in the problem are much shorter than the traveling distances 
between the interaction processes. In the models above, this can be easily 
achieved by violating Eq. (5). 

(2) All instable neutrinos have decayed before coherence is lost, which means 
that aL^"^ / E 3> 1. Since for Majoron decay a as well as L^"^ depends 
on the coupling constants, one can show that this condition together 




IQt^^E^ 



(4) 
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with Eq. (2) implies that £'^/Am^^ 3> 1 in the wave packet model above 
applied to Majoron decay. This is always true for relativistic neutrinos 
and small Am^'s, such as often assumed for active neutrinos. In addition, 
constraints such as Eq. (5) may have to be satisfied for the neutrino 
oscillations considered. 

Note that Sees. 4 and 6 are calculated for the first case. The last application, 
Sec. 7, corresponds to the second case. 



3 The formalism 

In this section, we will introduce the formalism used to calculate transition 
probabilities and fluxes. First, we will motivate the extension of the formalism 
in Ref. [12] by wave packet aspects and properties of point sources. Then, we 
will define the relevant operators, transition amplitudes, and fiuxes. Finally, 
we will give certain limiting cases in order to be able to describe realistic 
scenarios by simplified formulas and to show that the formalism reduces to 
the one in Ref. [12] in the coherent limit. 

3. 1 Motivation 

In Ref. [12], the combination of neutrino decay and neutrino oscillations was 
discussed for beams, but not for point sources. Baselines were assumed to be 
short enough such that the wave packets are still sufficiently overlapping at 
the position of decay and at the detector. Let us now generalize this with 
respect to the two above mentioned aspects. 

3.1.1 Coherence and wave packets 

In order to implement the loss of coherence by propagation over large dis- 
tances, L > L'^°^, we follow the mentioned wave packet treatment of neutrino 
oscillations in Refs. [22,23]. Therein it was noted that using wave packets 
leads to additional factors in the neutrino oscillation probabilities. The loss of 
coherence due to the spread and different mean velocities of sharply peaked 
wave packets is described by a factor exp[—{l/L'^^'^y] in the transition prob- 
abilities, where L'^'^ was defined in Eq. (3) and is the spatial width of 
the wave packet determined by the production P and detection D processes 
(cri)^ = {cr^y + {o-x Y- In addition, a factor enters the neutrino oscillation for- 
mulas, which is equal to unity if the constraint in Eq. (4) holds [23]. Violating 
this constraint formally also can be achieved by making the coherence length 



6 



very short. 

Taking into account an intermediate decay process X between production 
and detection, which corresponds to the detection process above, we expect a 
factor exp[— (/i/L^^*^'^)^] in the transition probabihties, where (cr^)^ = + 
{a^Y ~ (cr^)^ is the wave packet width for supernova neutrinos, li is the 
distance from production to decay, and the indices a and h refer to the mass 
eigenstates before decay. This factor describes the loss of coherence in the in 
state of a decay process, i.e., interference effects caused by fixed relative phases 
in an incoming superposition of mass eigenstates are destroyed. Furthermore, 
in decay new wave packets are produced. This leads to additional factors 
exp[— (/2/-^^c?''^)^] describing the coherence among the decay products c and 
d over the distance I2 from decay to detection. In this case, the wave packet 
width {(J'lY = {cr^y + (c^)^ is determined by the widths of the decay and 
detection processes. In Sec. 3.2, we will introduce an operator S, which has 
the suggested properties. 

3.1.2 Radially symmetric source fluxes and time dependence 

The second aspect we want to integrate into our operator framework is time- 
dependent, radially symmetric source fluxes. We define = dN^/dt to 
be the total source flux, i.e., the number of neutrinos of the flavor i>a emitted 
equally in any direction per time unit. Hence, ^^°^{t)dt = Na is the total 
number of neutrinos emitted by the source. In addition, let ^api^) be the "flux" 
at the detector, i.e., the number of neutrinos per time unit produced as flavor 
i^a and detectable as flavor up. Since this flux knows about the produced flavor, 
it could also be calculated by the transition probability as well as the number 
of Up arriving per time unit at the detector, i.e., Often the term/Zuxis used 
for the function However, in our notation — J2a^ai3- Furthermore, 
we will split ^a/si't) ^^^{t) and ^^^{t), describing neutrino fluxes with 
no intermediate decays and one intermediate decay such as the transition 
probabilities and P^^ in Rcf. [12]. These may either have to be added 
or not, depending on if the detector can distinguish undecayed particles from 
decay products (such as by energy resolution or spin) or not. 

Since decaying neutrinos, which are initially traveling into directions other 
than that of the detector, can produce neutrinos which are re-directed towards 
the detector, we average for the calculation of over all possible decay 
positions at which the secondary neutrinos may still arrive at the detector. 
Figure 1 defines the geometry of the problem, as well as the geometrical terms 
in the figure caption. It is safe to assume that the detector is small compared 
with the distance to the decay position a/D ^ I2. Thus, we can use radial 
symmetry, introduce spherical coordinates, and reduce the spatial averaging 
to / dcosx/2. For we only need to find the fraction of that can hit 
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L 



Fig. 1. The geometry of a beam from the supernova SN to the detector D. The 
direct baseUne has the length L, the baseUne with one intermediate decay has the 
length h + h, where the neutrino travels h before decay and I2 after decay. 

the detector by geometry, which is A{cosx)d/'^, i-^-, the cosine range under 
which the detector is seen by the source. 

Since a secondary neutrino produced by decay will arrive at the detector 
only with a certain probability due to kinematics, we introduce the function 

f]ij{L, li, I2, D) describing the fraction of decay products j of parent neutrinos 
i decaying at the position determined by /i, /2, and L, which will still arrive 
at the detector. It is given by 

77ij(L, li,l2,D) = r)ij{L, X,h,D) = ^ J 

ij p 

where T> denotes the surface of detection and D its area, which is to be inte- 
grated over, and Tij is the total decay rate. For the description of the kinemat- 
ics of mass eigenstates we will use the square root of this function, because 
amplitudes behave like square roots of particles [12]. We readily see from Fig. 1 
that l2{L,x,h) is determined by 

ll^ L"^ + ll-2Lhcosx, (7) 

i.e., we may, for example, choose li and I2 or x and li as sets of independent 
parameters. For Majoron models the redirection by decay is because of kine- 
matics limited by a maximum angle O^J^^^, where z/j is the parent neutrino and 
the decay product. This angle can be shown to be very small for Amf^- not 
larger than or by many orders of magnitude and relativistic neutrinos, 
which is reasonable at least for active neutrinos with no extremely hierarchical 
mass spectrum. Generalizing this to any considered decay model, using x and 
li as independent parameters, and assuming the detection area D to be small 
compared with the distance I2, i.e., \fD <^ I2, we determine the parameter 
ranges from geometry for small angles to be 



dV 



dcosO 



—d cos 9d(j), 



(6) 



xe[0,^j^J, (8) 

he [0,1-^], where ^ L (1 - . (9) 

V "max / 

Furthermore, we have to take care of time dispersion due to different traveling 
lengths. In general, oc - 1^'^ - 4 ), where ~ lk[l + ml/{2E^)] 

for relativistic mass eigenstates with mass nia, Vi is the mass eigenstate before 
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decay, and uj is the one after decay. Thus, in order to integrate the time 
dependence of the source flux, we incorporate it at the time when the peaks 
of the wave packets arriving at the detector were emitted at the supernova. 
Travehng times are related to mass eigenstates, which are, at least within 
the coherence length, to be summed over coherently. However, the supernova 
produces flavor eigenstates. Thus, it will become useful to introduce the notion 
of an amplitude flux v$^, describing the flux of states instead of particles. 



3.2 Operators, transition amplitudes, and fluxes 



We will now define the relevant operators, transition amplitudes, and fluxes 
based on the discussions in the last section as well as in Ref. [12]. The decay 
rate is defined as = ^ for i — > j decay, where is the rest frame lifetime 

for that decay channel. Time dilation by the factor 7j — Ei/rrii implies that 



T-iobserver rO ^,-1 ^ ^* '^^J 



The rate = a^j is the overall decay rate of the state i, which means 
that Bij = ^ is the branching ratio for i ^ j decay. The energy E^ of the 
mass eigenstate i will be approximated by the mean energy E of the produced 
flavor eigenstate if the energy corrections are of higher than first order. The 
operators will be defined in terms of effective creation and annihilation di 
operators operating on mass eigenstates [12]. 

Definition 1 (Disappecirance operator) is the transition operator which 
is also often called "decay operator". Effectively, it returns the amplitude for 
an undecayed state i remaining undecayed after traveling a distance I along 
the baseline L: 

V4l) = i:exp (-^] ala, with t« ^ (l + ^) /■ (11) 



Definition 2 (Appearance operator) is the differential transition op- 
erator, which destroys an in state and creates an out state in [ti,ti + dti] 
along the baseline L, i.e., a new state "appears": 




The probability density, which is the square of the amplitude, has to be inte- 
grated over li. Since i^'^ ~ (1 + ■^)l for a mass eigenstate i/j, the transition 
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amplitude, which is defined per time unit, needs to be transformed into an 
amphtude per length unit by the square root of this factor. Note that for 
a stationary problem the geometric function r] does not depend on traveling 
times. 

Definition 3 (Propagation operator) is the operator propagating a state 
Vk CL distance I along the baseline L: 

£\l) = exp {-lEkl) &%. (13) 



The index k will need to be summed over in the calculation of the flux, since 
time delays of different mass eigenstates will enter the (macroscopic) flux 
formula. 

Definition 4 (Decoherence operator) S^{l,T) is the operator describing 
the loss of wave packet coherence of a state by traveling a distance I along the 
baseline L: 

S^{l,T)^^N^''&q) 



{I - v,Tf 



Ma. 



r\2 



(14) 



It turns out that this operator has the required properties leading to the 
factors postulated in Ref. [23] as well as in Sec. 3.1.1 in order to describe the 
loss of coherence for long basehnes. The variable T describes the wave packet 
distribution in time. It will be integrated over in the transition probability, 
since it is not measurable by the target process. In addition, o"^ is the wave 
packet width of the composition / of two processes (production, decay, or 
detection). Furthermore, iV^'* is a normalization factor with 



/ 



N^''N^'^ exp 



(/ - v.TY + (/ - v.Tf 



dT = exp 



' coh,7 



(15) 



Here the factor exp{— (/ — faT)^/[4(o"^)^]} comes from the expansion of a 
sharply peaked wave packet, which can be written as ^a(-^i T) = exp[—iEaT + 
ipaL — (L — faT)^/(4(T^)]. In the definition of the propagation operator we 
ignored the factor exp(ipaL), since it would only give rise to common phase 
factors as well as an additional, already mentioned factor, which is negligible 
for (7^ < L°^^ 

Definition 5 (Calculation of transition amplitudes) The transition am- 
plitude A{vo, — > vp)'^ — A^p by exactly n intermediate decays is for n — 



A%^{vp\S\L)V_{L)S^^mT)W) 



(16) 
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and for n — 1 

(17) 

The time variables T and U are to he integrated over in the transition prob- 
abilities, i.e., after taking the squares of the amplitudes. The indices i and j, 
denoting the intermediate traveling mass eigenstates, are to he summed over 
in the calculation of the fiux. Here (cr^^)^ = (cr^)^ + (c^)^, where the label 
P refers to the production process, D to the detection process, and X to the 
decay process. 

Note that the order of T>_, S, and S is arbitrary, since it can be shown that 
these operators all commute. 

Definition 6 (Calculation of fiuxes) The supernova neutrino fluxes ^^^^ (t) 
(no intermediate decays) and^^^{t) (one intermediate decay) at the detector 
D are calculated as 



0,1 



i 

-I /" 



COS^max 'l"*^ 



» 3 



dh 



dcosx 



(18) 
(19) 



with r« ~ (l + ^) L, i« ~ (l + ^) I,, 0^^ = max^^^, and = 
maxZ™^. 

In we must integrate over the maximal possible range determined by the 
maximum of all O'^J^^ and l^^f^, in order to cover the appropriate spatial region. 
If (x, li) is out of range, which may occur in some integrands, the function rjij 
in X>+ will evaluate to zero. 



3.3 Limiting cases 



We now discuss some limiting cases. Most of them involve coherence limits 
and can also be combined with the limiting case in Sec. 3.3.3, i.e., neglecting 
different travehng path lengths. 
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3.3.1 Coherent wave packets 



Let us assume that the wave packets are at all times coherent, ie., L <C L'^°^ 
for all coherence lengths involved in the problem. Hence, the decoherence 
operators S^{l,T) will only give rise to factors of unity after the integration 
over time T, and hence, can be neglected. In addition, this implies that Ik — tk, 
i.e., we do not have to take into account differences between traveling times 
and propagation distances of different mass eigenstates. Finally, we assume 
that 91^^^ ~ 6'max and l^f^ ~ Z™*^"^ to be independent of the indices i and j, 
i.e., within the coherence length the mass eigenstates take approximately the 
same traveling paths. Therefore, we can pull the amplitude fluxes out of the 
summations over the intermediate traveling mass eigenstates, and redefine the 
propagation operators by absorbing the summations. In addition, we introduce 
probabilities instead of amplitudes. We then obtain the following expressions 
for the operators: 



"-^ aja,, (20) 



V+(L,h,l2,D)^Y.T.\Rr ^Vij{L,h,l2,D)a%, (21) 
£{l)^^e^p{-iEil)alai. (22) 



For the transition probabilities we have 



P'a^ = \{Mm'r>-iL)K)\\ (23) 
= I {up\£{k)V4k)V+{L, h, k, D)£{h)V4h)\u^)\' , (24) 

and for the fluxes we find 



^^Jit)-^l'\t-L)^^^^P^,, (25) 
/ (26) 

COS^max h=0 

Therefore, the expressions found in this limit are very similar to the ones in 
Ref. [12], but adopted to point sources. 



12 



3.3.2 Incoherent wave packets 

This limit corresponds to the first case mentioned in Sec. 2, i.e., loss of co- 
herence between any two processes in the problem. Thus, we assume that 
L, h, I2 ^ L'^"^ for all coherence lengths in the problem. Hence, the operators 
S\l,T) will give rise to factors exp[— (//L™'^''^)^] — > 8ij after the integration 
over the time T. It can be shown by application of the operators that this cor- 
responds to incoherent summation over the intermediate traveling states, i.e., 
iXlij/ijP = Hij\fij\'^ in this limit. Therefore, we can square the amplitude 
fluxes, put the summations in front of the integrations in and contract 

the integration limits back to the appropriate regions. We then may define 
transition probabilities 

i^a;^=l(^/.k.)l'l(^d^^-Wk.)n(^.k.)l% (27) 

^=i(^/^i^.)n(^.i^^-(^2)p+(i.,/i,/2)2^-(Mk.)n(^.k.)i', (28) 

and fiuxes 



*S°(^) = E ^Tit - T») pO. A(cosx)D ^ (29) 

pnax 

^^/(t)=EE / ] ^':\t-tf-tf)'^-^di,^-^ (30) 

with T(*) ~ (^1 + L and t'i^ ^ (l + h- Thus, S splits up the formu- 
las for the transition amplitudes such that neutrino oscillations vanish. 



3. 3. 3 Neglecting time delays by different traveling paths 

Time delays due to different path lengths can be neglected when they are small 
compared with time delays by different masses. The fraction of the initial flux 
arriving at the detector is D/AttL"^, since the total flux through any sphere 
around the point source is equal by symmetry. This is, in principle, not changed 
for decays from one mass eigenstate into another, since the overall number of 
neutrinos arriving at the detector is unchanged.^ Nevertheless, the relative 
arrival times may differ for different path lengths. The time delay due to 
different path lengths Ati^ij for mass eigenstates z/j and uj, propagating before 
and after decay, respectively, can for small 91^^^ be approximated by Ati^ij < 

^ Note that we only assume small changes in direction by decay, which means that 
the detector geometry does not affect this conservation law. 
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^Vf^J ~ L{9]ig^y/2 from geometry. For this small time interval, delays by 
different masses can be ignored as second order effects. In order to have the 
limiting case discussed here, it needs to be much smaller than the time delay 
by different masses At2,ij = Am^jL/{2E^). Therefore, for 6'^^ < JAmfj/E 
effects of different path lengths can be neglected. In fact, it can be shown 
that this is equivalent to Amfj <^ rrij in the case of Majoron decay. However, 
the time delay caused by different path lengths may still be measurable if the 
absolute value of Ati is longer than the time resolution of the detector. 

Since we are ignoring traveling times due to different paths and since the total 
number of arriving particles is not changed by decay, we can implement this 
limiting case in the geometrical function rj, describing the fraction of secondary 
neutrinos hitting the detector, as 

?7(L, X, k, D) = ^ 4(5(1 - cos x)- (31) 

Thus, the secondary neutrinos are peaked in the forward direction and for 
any x < \ the integral /J rj d cos x/ 2 = D/{A'kL'^) gives the required fraction 
independent of the indices i and j. Note that this definition of rj corresponds to 
a forward peaked differential decay rate dVij/dcosO oc 5{1 — cosO) in Eq. (6). 
However, it is not exactly identical due to re-direction effects. Nevertheless, 
a distribution of the differential decay rate, which is sharply peaked into the 
forward direction, gives an effective 9li^^ ^ Oli^^. Then, for the case of dli^^ <^ 

' O illtiX ill tlx ■ lllciX 

Amfj/ we have this limit again. We will introduce in Sec. 5 an r;-function 
which corresponds to this case. 

In order to incorporate the new r^-function in this limiting case, we only need 
to re-define one operator: 



^+ = EEJ^\1 + :S^«K (32) 




In addition, we need to take into account that I'^ — L — li'm. the evaluation of 
the ^-integration over the ^-distribution in the 77-function. Finally, we obtain 
for the fluxes for /i = / 



0,j 



D 



47rL2' 



dli 



D 



47rL2 



(33) 
(34) 



with r» ^ (1 + ^) L, ^ (1 + S) I: and ^ (l + ^){L-l). 
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Simultaneous application of the incoherent limit yields in addition formulas 
adjusted to the last limiting case in a straightforward way. 



3.3.4 Decay before loss of coherence 

In this limit, we treat the second case mentioned in Sec. 2, i.e., decay rates 
large enough such that all of the neutrinos decay before coherence is lost. If 
we assume that L > L""^, the wave packets of undecayed states will loose 
coherence before detection, i.e., can be calculated as in Sec. 3.3.2. For 
very early decays, i.e., a ^ E/L^°^, the distance I2 is much longer than the 
corresponding coherence length, whereas li is shorter. In this case, all initial 
neutrinos have decayed before li approaches the coherence length, which also 
means that ~ 0. Furthermore, since L » L^°^, the decay products will 
loose coherence before detection. In order to understand what may happen in 
this limiting case, we assume interference between different decay channels, 
i.e., simultaneous couplings to the decay products, such as by Majoron cou- 
pling constants > and gjk > for different indices i, j, and k. Then, 
the relative phase of the states in the incoming superposition depends on the 
decay position. In addition, the arrival time depends on the decay position, 
because the mass eigenstates travel with different velocities before and af- 
ter decay. We expect to observe the most interesting interference effects for 
Am^ ~ a. For fast oscillations, i.e., Arn^ S> a, interference effects will be 
washed out by averaging over the different decay positions. For slow oscilla- 
tions, i.e., Am^ <S a, all particles will still have the initial relative phase at 
the decay position, eliminating any oscillating effect. 

In order to calculate the transition flux we combine the formulas in 

Sees. 3.3.1 (before decay) and 3.3.2 (after decay). We define 



^'^[^^2E~^)^' ^^^^ ^'^]7E< (35) 

where the mass square average is to be taken over the N mass eigenstates 
oscillating before decay. Since this mean traveling time will only enter in the 
total source flux and li is assumed to be quite short, this is certainly a reason- 
able approximation for the source flux not changing too much on timescales of 
the order Am^/i/(2£'^). In addition, we assume O^i^^ ^ 6^^^^ and ['{"^f^ ^ l^j"" to 
depend only on the decay products (the incoming wave packets are coherent). 
Eventually, this yields 
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dli 



incoherent 

/max 
1 l,j 



coherent 



^ cos ei,,, '1=0 



dh 2 ' 

with the operators from the coherent part as given in Sec. 3.3.1. 



(36) 



(37) 



4 Incoherent mass dispersion 



In this section, we will demonstrate for the limiting case of incoherent propa- 
gation how mass dispersion can be caused by decay even if we ignore different 
traveling path lengths.[^ This is done by simultaneous application of the lim- 
its in Sees. 3.3.2 and 3.3.3. We will show that it is possible to observe a 
time dispersion by decay even for a pulsed source flux. Thus, we assume that 
$a°*(t) = Na S(t + L) which describes a neutrino pulsed produced at t = — L, 
where is the total number of neutrinos of flavor emitted by the super- 
nova. In this case, it will be possible to detect massless neutrinos at t = and 
massive neutrinos delayed by At*^*) = Decay leads to a time dispersion, 
since a mass eigenstates travels with a different velocity before and after decay. 
Its arrival time depends therefore on the decay position. 

Combing the limiting cases in Sees. 3.3.2 and 3.3.3 as well as applying the 
transition probabilities and operators, we obtain for 



2 



The 5-distribution in this formula, describing the neutrino pulse, implies that 

Am,^. rn?- , ^ 

t = ^/ + -^L 39 

2^2 2^2 ^ ^ 



^ A similar effect was discussed in Ref. [28] for the arrival times of photons as decay 
products of radiative neutrino decay. 
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or, for rrii ^ rrij, 



The condition < / < L, for which the ^-distribution evaluates to a non-zero 
value by the integration, is then equivalent to 

§L < < ^,L. (41, 

Thus, the signal arrives at the detector between the signals of the undecayed 

heavy mass eigenstate and the light one produced by the intermediate decay. 
Note that the allowed time interval depends on the indices i and j. In ad- 
dition, we need to take into account that 5 (/(/)) = 6{l — a)/|/'(a)|/(a)=o = 
2E'^ / AmfjS{l — a) with a equal to I as given in Eq. (40). Finally, integrating 
over I leads to 



2 ,2 ^Eaij ( mj 



X exp 



1 + 



2E^ 



- a,: 



1 + ^ 
2E^ 



2E 
Ami 



(42) 



In the limit of stable decay products, i.e., aj —>■ 0, and by ignoring corrections 
of the order m'^/E'^ to the signal height, this reduces to 



2Ea,. 



Amfj 



2^ exp 



-ai 



2E 



t - 



2^2 



(43) 



Thus, the flux pulse of the source is smeared OTit to an exponentially dropping 
signal at the detector. Here the allowed range for tij depends on the mass 
eigenstates i and j and is given in Eq. (41). One can immediately see that the 
smearing is determined by the coefficient a^E / Amf^ in the exponential. We 
obtain maximal smearing for small aiE/Amjj, i.e., small a^'s or large Amf^s. 
A numerical analysis for typical values of Amjj, E, and L shows that the 
aij's, and the Majoron coupling constants gi/s, respectively, can be far below 
any currently assumed upper limit for observing exponential dropping of this 
function (for constraints on the Qi/s, see Refs. [29-32]). However, if ai is too 
small, the factor aij will suppress the term completely. 

We illustrate the effect for the following three-neutrino scenario: Maximal 
mixing, i.e., lua) = 75^1) + 75^2) + 75^3) (6*12 = On = 6*23 = 45° and 
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5cp = 0), decay of v-^ into z/2 or Vi only, i.e., 032 = 031 = 10 ^^GcVkm 
E = lOMeV, L = 10^2 m ~ 32kpc, = 9 ■ (AtiL^)/D, = 4cV, ms = 
2 eV, and mi = 1 cV, which means that Am|2 = 12 eV^ and Am|i = 3 eV^. We 
can use Eqs. (18) and (43) to evaluate and where the sums are split 
up in order to see the different signals from different mass eigenstates. Figure 2 
shows the separated signals of the decay products as well as the ones from 
the undecayed particles ^^f- been observed [33] that the SN1987A data 
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Fig. 2. The (separated) signals of the decay products and those of the unde- 
cayed particles ^ai^ for the scenario constructed in the text. The (^-distributions in 
^af are plotted as Gaussian signals, where the leftmost refers to the lightest and 
the rightmost to the heaviest mass eigenstate. The vertical line at t = indicates 
the reference time when massless neutrinos would arrive. 

can be reasonably fit to a decaying exponential with time constant r ~ 3s. 
Figure 2 indicates that we can easily find parameter sets in order to have 
an effect with a similar time dependence. However, for the currently assumed 
parameter values for active neutrinos, as shown in Fig. 3, the timescale of 
this effect would be rather short. In addition, the two different decay channels 
could, mainly due to the small Am|i, not be resolved anymore. Thus, in 
order to be able to identify all relevant individual effects in the plots, we will 
further on choose parameter values which amplify these effects more clearly, 
but which are not the values currently favored by data. There are basically two 
arguments supporting this point of view. First, the best-fit parameters may 
change when making global fits with neutrino decay included, and therefore, 
the mass squared differences may become larger. Second, more recently, other 
models being relevant for this sort of discussion have been introduced, such as 
in Ref. [34] , in which the decay of a light sterile neutrino into active neutrinos 
was proposed. In this case, the relevant mass squared differences could again 
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Fig. 3. The signals of the decay products and those of the undecayed par- 

ticles ^ai^ for the scenario with the parameter values Am|2 = 2.5 • 10~^eV^, 
Am|^ = 3.65 ■ 10~^ eV^, mi = 10~^ eV (absolute neutrino mass value needed here), 
^12 = ^23 = 45° (bimaximal mixing), ^13 = 5°, and Scp = 0. The other parameter 
values are the same as for Fig. 2. Note that a can be equal to /i or r. For a = e the 
result will be different, since the mixing is not trimaximal anymore as in Fig. 2. Fur- 
thermore, note that the scale of the fluxes is not applicable to the 5-fluxes illustrated 
as peaked signals, but the heights of these peaks reflect their relative values. 

be rather large. Such light sterile neutrinos could be produced by oscillated 
supernova neutrinos close to or within the supernova. 



5 Dynamics of Majoron decay 



Let us now give an approximation for the function r]ij{L,x,li, D) defined in 
Eq. (6) for Majoron decay of relativistic neutrinos. It describes the fraction 
of decay products arriving at the detector for decay at (x, /i) (c/.. Fig. 1). 
Since only a small fraction of neutrinos will decay close to the detector, we 
can assume that 



where 6^J^^ is given by [12] 



L > » /D/e«„, (44) 
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Here the index i refers to the parent neutrino and j to the decay product. 
Hence, we can treat the detector as a point target and approximate Eq. (6) 

by 



dVij 



dcos9 



(cos 9)d 



1 

2^* 



1 

J- ij 



dr 



dcos9 



{cos6)d 



D 
2^ 



(46) 



with I2 given by Eq. (7). From 



Ei — Ej 



V>i? + \v>i? - 2|Pi||Pj| cos6l, 



(47) 



which determines the kinematics of the process, we find for small 9 two energies 
Ej corresponding to one angle 9 



Ef{9) = 



Eiil+x^ 



^ + X' 
j 



1 ± 



m. 



'\\X 



2\2 



(48) 



Equation (46) can be written as 
1 D 



dEj 



E1 



dcos9 



+ 



dTii 



ET 



dEi 



dcos9 



E- 



(49) 



The differential and total decay rates in this equation are for (pseudoscalar 
or scalar) Yukawa couplings in the Lagrangian for decay into neutrinos or an- 
tineutrinos given in Refs. [12,27], respectively, for the Lagrangians introduced 
there. Prom Eq. (47) we read off 



dE, 



dcos9 



E^ 



iPillPjl 

\Ei — Ei 



E, 



- 1 



(50) 



Since we use x h as independent parameters, we need to express 9 in 
terms of x- To first approximation for small angles, which is reasonable for a 
small 9'^J^^, we find from geometry 



9 



^1 



L-h 



(51) 



Now 77 can be evaluated numerically and is plotted in Figs. 4, 5, and 6 for 
(pseudoscalar or scalar) Yukawa couplings in the Lagrangian and decay into 

neutrinos and antineutrinos. Note that the parameter ranges determined in 
Eqs. (8) and (9) imply that r] evaluates to zero, i.e., nothing arrives at the 
detector if the parameter set lies within the upper-right triangle over the di- 
agonal. In this case, the secondary neutrino cannot be re-directed towards the 
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Fig. 4. Density plot with a linear shading scale of the function f]ij{Xjh) for 
pseudoscalar Yukawa couplings in the Lagrangian (decay into neutrinos). Here 
D = 10000 m^, L = 10^2 m ~ 32kpc, Ei = E = lOMeV, = 2eV, and mj = 1 eV. 
For the geometry of the problem, see Fig. 1. 



detector anymore for kinematical reasons. The plot for antineutrinos shows, in 
comparison to the neutrino plots, that the direct path (x = 0) is suppressed 
for small li, because of a spin flip. For Yukawa scalar couplings 77 turns out 
to be quite independent of the angle x- cases, we can expect heavy 

suppression for too large x- 

The magnitude of 7] is mainly determined by the D/Z^-dependence of Eq. (46) 
(antineutrinos) or by the factor dEj/dcos9 in Eq. (49) (neutrinos). In our ex- 
ample, D/ll — 0{D/L'^) — 10~^° for large I2. One can show that dEj/dcos9 
diverges for Ej = E^ at ^ = 0, which evaluates there to E'^ = Ei. Thus, 
in the limit L > I2 ^ V^/^maxi the ^^-function is difficult to evaluate nu- 
merically in the transition probabilities. Nevertheless, since the integral over 
the differential cross section over the whole parameter range has to give 
Fjj, the mean differential decay rate is of the order ^ij/{Oll^^)'^, and thus, 
fj — O (£)/(Z2^J^ax)^) ^ -^/^2- Hence, for neutrinos the main contribution to 
the differential decay rate comes from the divergent, forward peaking part, 
and we may, in this case, neglect time delays by different path lengths as it 
was done in Sec. 4. For antineutrinos the forward direction is suppressed by 
the spin flip and becomes finite, which means that we may approximate the 
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Fig. 5. Density plot with a linear shading scale of the function i]ij{Xih) for scalar 
Yukawa couplings in the Lagrangian (decay into neutrinos). Here D = 10000 m^, 
L = 10^2 m ~ 32kpc, Ei = E = lOMeV, = 2eV, and ruj = leV. For the 
geometry of the problem, see Fig. 1. 

differential cross section by its mean value in order to obtain from Eq. (46) 
fjiiiL, X, h, D) = — ( -^^] —An = . (52) 

r ■ ■ 

Here O^l^^ < ^^^lax) where is an effective maximum angle. For neutrinos, Figs. 4 
and 5 indicate that 6*^^^ <^ (^max could also be used as an approximation. For 
antineutrinos. Fig. 6 is dominated by the l/Zl-dependence, which means that 
we may use the approximation ^^^^ = ^^Aax- Note that for ^^g^ < ^J^^^ we 
also have to change 9^^^^ ^ 9^^ and l^^ ^ l^^ = L(l - x/OHJ in the 
integration limits in Eq. (19). 



6 Fluxes for incoherent propagation 

We will now provide the most general expressions for the fluxes in the limit 
of incoherent wave packet propagation introduced in Sec. 3.3.2, where we also 
take into account time delays by different path lengths. Again we assume a 
neutrino pulse produced at the supernova, i.e., $a*(t) = NaS{t + L), such 
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Fig. 6. Density plot witli a linear shading scale of the function r]ij{x^h) for pseu- 
doscalar or scalar Yukawa couplings in the Lagrangian (decay into antineutrinos) . 
Here D = 10000 m^, L = lO^^ m ~ 32kpc, Ei = E = lOMeV, mi = 2eV, and 
rrij = 1 eV. For the geometry of the problem, see Fig. 1. 

that massless neutrinos would arrive at the detector at time t = 0. For other 
sources the results for this 5-flux can be used to describe a time-dependent 
source. We neglect corrections to the signal heights coming from / ^ t, since 
these are of the order vr? jE"^ . However, we are interested in time dispersion 
induced by non-zero masses. Noticing that 27rA(cosx)£)-^^ = -D for the direct 
baseline, we obtain from Eqs. (18) and (27) 

^ll^J - ^l^a.nC/,.re-^^ (53) 

Prom Eqs. (19) and (28) we find 



^max 

V/ coses.. '1=0 ^ ^ ^ ^ ^ ^ 

X |C/.,HC/;3,|V^'^^77,,(L,x,Zi,i^)e-^'MZi^^. (54) 
hi z 

^ V ' 
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Integrating over /i, wc must observe that the hmits of the Zi-integration depend 
on X s-iid we have to take into account that the (5-distribution evaluates to a 
non-zero value only for < /i < = L(l — This can often be 

done by adjusting the integration limits of the x-integration, but here it will 
lead to quite complicated expressions for li{t, x) and x(^, ^i), respectively. We 
therefore introduce the function 

f 1 for < hit, x) < Ifff = l(1- 
I otherwise 

to describe the region where the integrand contributes. Furthermore, we have 
to take into account the transformation of the (^-distribution 6 (/(/)) — ^{l — 
a)/|/'(a)|/(a)=o with a — li{t, x) being the solution of /(a) = (see below). In 
this case, we can approximate |/'(a)|/(a)=o by 

|/'(a)U,,.„.^ + ^j^^f (56) 

Thus, for small x we arrive at 



^+3 



Cmax c (+ \ 

X / e-i("^'^+"^'^)r^,,(L,x,ii,i^)^^X^X (57) 

with 

^1 = ^ ^ (58) 

coming from the 5-distribution. In addition, Eq. (7) implies that l\ — L'^ + lf — 
2Lli cosx- This leads together with Eq. (58) to a quadratic equation in li{t, x) 
or l2{t,x)- Analysis of li{t,x) shows that we obtain a unique, quite lengthy 
solution for li > 0, which we will not present here. It can also be shown 
that t(x, h) grows monotonously with growing x or /i, which implies that the 
earhest arrival time is Lm^/2£'^. The latest arrival time can be obtained from 
the maximum of t{x, IT^ix))^ which is again a quite lengthy expression. Since 
we neglect corrections of the signal height of the order w? /E'^ and second 
order corrections, we can approximate Eq. (58) in the exponential of Eq. (57) 
by 

k-t-h + L. (59) 
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We finally obtain 



^+3 



xr7.,(L,x,/i(t,x),^)^|^Xt^X- (60) 



For an arbitrary time-dependent source = ^af^ + -f') with fit + 

L)dt — 1 we can calculate the flux at the detector from the above expres- 
sions for and by 

oo oo 

^ap{t)= J ^Sp{t-t')f{t' + L)dt' = — J ^^,{t-t')^T{t' + L)dt'. (61) 

— oo — oo 

Similarly, one can fold these fluxes with energy dependencies and detector 
properties. 

Let us now illustrate the effects of decay in Eq. (60) by an example similar 
to the one at the end of Sec. 4. Wc are using the same parameters except 
from asi = and 0:32 = cao = c 10^^^ GeVkm^^, where c = const., since 
one decay channel is sufficient for showing the effects. For rj we use the ap- 
proximation fj in Eq. (52), which also means that the upper integration limit 
in Eq. (60) is to be replaced by 9l^^^ = 9^^^ < 9f^^^ = 6'max, as well as 9l^^^ 
in Eq. (55), which represents the integration limits for li. It turns out that 
the numerical evaluation of Eq. (60) is quite sensitive to the approximations 
as well as to the parameter sets, which means that the solutions only can be 
used to demonstrate the qualitative behavior. From Eq. (55) we know that 
the function ^{t, x) can only take the values or 1, giving the areas where the 
integrand in Eq. (60) is defined. Figures 7 and 8 show these areas times the 
traveling path length before decay li{t,x) for i9max = 6'max and 9^^^ = 9^^J2. 
Note that the x-scales are different in these two plots. The white regions 
indicate ^ = 0, i.e., the integrand in Eq. (60) does not contribute. One can see 
that the contours of equal li become quite independent of x small ^^max, ^-c, 
the arrival time t becomes dominated by different velocities due to different 
masses (c/.. Fig. 7). For large 6'max we can see the effects of the different path 
lengths, curving the contours towards large t, ie., late arrival, for large x (c/.. 
Fig. 8). Thus, in comparison to the example in Sec. 4, we expect enhance- 
ments for late time arrivals and suppressions for early time arrivals, because 
of the delays due to different traveling paths. For the case of Majoron decay 
in Sec. 5 we demonstrated that for neutrinos the function ?7(x, h) is favoring 
small X) which means that the approximation ^max <^ ^max should be good. 
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Fig. 7. The linearly scaled contour lines of the function ^{t, x) h for ^max = d„ 
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Fig. 8. The linearly scaled contour lines of the function ^{t, x) h for ^max = ^max/2. 

In this case, the problem reduces to the one in Sec. 4. For antineutrinos, the 
approximation Oiaax — ^max is better. 

Figure 9 shows the quahtative behavior of for different O^ax- The verti- 



26 




t [s] 



Fig. 9. The flux of decay products $Q,e for the decay rate a = ao = E/L and 
different values of ^max- 



cal linesP^ indicate the arrival times of the light and heavy mass eigenstates 
on the direct paths, described by the (5-distributions in whereas mass- 

less particles would arrive at t = 0. One can see that the smaller ^max, the 
more exponential the time dependence becomes, such as in Sec. 4 for forward 
traveling only. For large ^max late time arrivals are preferred, corresponding 
to particles delayed by longer path lengths. In principle, some particles could 
even arrive after the heavy mass eigenstate traveling on the direct path, but 
these had to decay quite late in order to keep the slow velocity of the heavy 
mass eigenstate as long as possible. Geometry {6 < ^max) implies, however, 
that we have the shortest overall traveling paths for li close to or L. This 
means that this is only possible for the case where path length effects domi- 
nate the mass dispersion, i.e., fi'J^^x ^ y'Amfj/E', which we do not consider 
in this example. 

In Fig. 10, the effect of different decay rates on is illustrated. Note that ao 
is chosen such that aoL/E = 1. For larger decay rates we see a behavior closer 
to the one in Sec. 4, i.e., exponential dropping. In this case, most particles 
will decay early along the traveling path. Geometry implies in addition that 
the overall path length is close to its minimum for li ~ 0. Therefore, for large 
decay rates path length effects can again be ignored. We conclude that different 
traveling path lengths only have to be taken into account in computing the 



^ The functions are not plotted in the whole possible range, because close to the 
limits the numerical evaluation becomes quite unstable, since we need to integrate 
an almost divergent function over an infinitesimally small interval. 
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Fig. 10. The flux of decay products for ^max = ^max and different values of 
decay rates a as multiples of ckq = E/L. 

supernova neutrino signal for large enough ^max and small enough decay rates. 
However, since is directly proportional to the decay rate, too small decay 
rates will make the flux $£;^ vanish. 



7 Early coherent decays 



As already mentioned above in Sec. 3.3.4, one can sec some peculiarities for 
incoherent propagation over the entire baseline L, i.e., L ^ but for 

decay before loss of coherence, a ^ E/L^°^\ Let us again assume a flux pulse 
^tot _ 7Vo,5(t + L). Neglecting corrections to the signal height of the order 
m^L/(2_E^), as well as path length dependencies for simphcity, i.e., combing 
the results of Sees. 3.3.3 and 3.3.4 in a straightforward way, we obtain from 
Eqs. (34-37) 
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where refers to the mean mass square of the heavy mass eigenstates before 
decay, as introduced in Sec. 3.3.4. Defining K'^j^i = U*iUi3jUakU^i and assum- 
ing the secondary neutrinos are stable, i.e., aj = 0, as well as '^K = 0, we 
find after some algebra 



k j 



47rL2- 



X exp 



ijkj 



t - 



cos 



t - 



rri; 



(64) 



Here m~}. = {ml + mi)/ 2 is evaluated for the heavy mass eigenstates before 
decay, and the condition m|L/(2£'^) < t <^ m}f.L / {2E'^) , coming from the 
integration limits, has to be satisfied. This is a very interesting result, since 
we may see a time-dependent oscillation if two ccy's are non-zero, i.e., if we 
have two decay channels with only one neutrino decay product. The time t 
indirectly measures the decay position via the fiight time of the decay product. 
Thus, it also measures the relative phase of the incoming states at this position, 
which induces the oscillations. However, since all a's are assumed to be quite 
large, a ^ E/L'^"'^, the exponential indicates that the dominant contribution 
comes from t ~ m^L/{2E'^). Therefore, only a small fraction of the possible 

time interval At 



ik 



mJ)L/{2E'^) will be covered by the signal. The 
oscillating efi'ect will hence not be observable in most cases. Thus, we have to 
integrate the observed signal over time, leading to 



i k j 



\OLi 



akf + (Am?,) 



In comparison to entirely incoherent propagation without interference of decay 
channels, which corresponds to bik in the summation, we observe the following 
two peculiarities: 

• Interference of decay channels aij ^ 0, akj 7^ 0, and i ^ k creates addi- 
tional terms in the summation, which, in principle, enhance the number of 
detectable neutrinos. 

• The mentioned interference terms are reduced with increasing Amfj^ due to 
phase shifts in early decays. This reduction can be neglected for Amf, <C 
Oi + ak, because the decay length is in such a scenario much shorter than the 
oscillation length and the neutrinos will have decayed before any oscillations. 
For Am^j. ai + ak the interference terms will vanish by averaging over all 
possible decay positions and the corresponding phases. 
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8 Summciry and conclusions 



In this paper, we have combined neutrino decay and neutrinos oscillations 
for radially symmetric sources, such as supernovas. We have calculated time- 
dependent fluxes at the detector by taking into account decoherencc due to 
the long baselines, interference effects within the coherence lengths, time de- 
pendence of the source, different flight times of different mass eigenstates, and 
different travehng paths for neutrinos re-directed in decay. We observed two 
interesting arrival time dispersion effects, which can in some cases be of the 
same order of magnitude: 

(1) The time of arrival depends on the decay position, since for visible neu- 
trino decay products neutrinos travel with different velocities before and 
after decay. The decay positions are distributed in an exponential man- 
ner, which means that we will also observe an exponential decrease of the 
detected flux in time. 

(2) The time of flight depends on the path. For radially symmetric sources, 
neutrinos with different traveling paths arrive at the detector, since de- 
cay may change the direction of the secondary neutrinos. This effect sup- 
presses early time arrivals and favors late time arrivals, working in the 
opposite direction to (1). 

We have demonstrated that the flrst effect can mimic an exponentially de- 
creasing flux at the detector, such that it can be fit to SN1987A data, for very 
small decay rates. This shows that supernova properties, which are inferred 
from a description without neutrino decay, can be altered or even completely 
changed. Moreover, we have shown for Majoron decay that the second effect 
has only to be taken into account for antineutrinos and not too large decay 
rates aL/ E ~ 1, if the masses of the participating neutrinos are of similar or- 
der of magnitude. For the case of Am^^ 3> for active neutrinos, the second 
effect would be much larger than the one discussed in this paper. 

Coherence is not only related to the production process, but also to the detec- 
tion or decay process. Thus, for small coupling constants in the Lagrangians 
coherence lengths can in some cases be quite large. In some sense, neutrino 
decay may act as a "coherence lens" by making the wave packets collapse. 
However, finite lifetime constraints, such as Eq. (5) for Majoron decay, have 
to be satisfied in order not to wash out interference effects. There are also 
some counter-intuitive peculiarities, which may be observed for very large de- 
cay rates. Even if only incoherent mass eigenstates (or at least one) arrive 
at the detector, there may be interference effects modifying the event rates. 
This is because neutrinos for very early decays may decay as long as they 
are still coherently propagating. For extreme choices of the masses and decay 
rates, as well as simultaneous coupling to the decay product by two differ- 
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ent mass eigenstates in the Lagrangian {e.g., for Majoron coupling constants 
g^i ~ (721 7^ 0). one can, in principle, have a time-dependent oscillation of the 
signal at the detector. In this case, the sensitivity-dependent neutrino oscilla- 
tion is transferred into a time dependence of the signal by the time of flight 
concept of the different mass eigenstates. Note that this may also happen even 
if only one mass eigenstate is stable, which flnally arrives at the detector. 

Since extremely small coupling constants in the Lagrangian can cause neu- 
trino decay over typical supernova distances, neutrino decay should always be 
taken into account in the calculation of transition probabilities and fluxes. The 
observations of supernova neutrinos from SN1987A indicate so far only that 
there is at least one stable mass eigenstate. We have shown that a number of 
effects involving neutrino decay into different mass eigenstates may alter the 
event rates at the detector even for only one arriving mass eigenstate. Thus, 
supernova neutrinos are an excellent probe to test neutrino decay. However, 
similar effects can be induced by details of the supernova explosion, which 
means that a separation of decay effects and supernova details may be diffi- 
cult. 
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